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Presented herein Is a parametric study of the defining equations which 
govern the steady state operational characteristics of the Grumman Monogroove 
Dual Passage Heat Pipe. These defining equations are combined to develop a 
mathematical model which describes and predicts the operational and performance 
capabilities of a specific heat pipe,given the necessary physical characteristics 
and working fluid. 

Included is a brief review of the current literature , a discussion of the 
governing equations, and a description of both the mathematical and computer 
model. Rnal results of preliminary test runs of the model are presented and 
compared with experimental tests performed by Grumman on actual prototypes. 
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INTRODUCTION 


Presently, the primary means for rejecting heat from orbiting spacecraft 
is through a space radiator system composed of a series of fluid loops. These 
loops circulate fluid through the radiator panels which in turn reject heat 
to the space environment. Because the current system uses a mechanically 
pumped coolant circuit to transfer heat throughout the radiating surface, it 
results in a system whose long mission reliability is low and one vulnerable 
to complete failure due to penetration by a single meteoroid. Reliability can 
be increased through the use of redundant plumbing, pumping, and valving hard- 
ware, resulting in a large increase in total system weight. Hence .there is a 
need for significant technical improvements in the development of a long life 
heat-rejection system which is suitable f ' long term .high power missions and 
can be constructed and deployed on orbit. 

One solution to this problem is the development of a large modular radiator 
system that can be assembled during orbit from a number of standard components. 
This space-constructable radiator system would fulfill the needs and demands of 
large long-lived heat rejection systems and would allow systems to be built up 
to any desired heat load capacity. 

The key component of this concept as it is presently conceived is an inno- 
vative. high-capacity, dual passage heat pipe designed by Grumman Aerospace 
Corporation (1). This heat pipe with radiator fins attached, would be "plugged 
in" to contact heat exchangers providing heat removal from a centralized heat 
transport looo. This type of system would be insensitive to complete failure 
due to micrometeoroid puncture, with the puncture of any single heat pipe re- 
sulting in only the loss of that module's 2-kilowatt capacity. The damaged seg- 
ments could be removed by the Orbiter and replaced or repaired as necessary. 

The basic design of this improved high-performance dual passage, heat-pipe 
consists of two large axial channels, one for vapor flow and another for liquid 
flow (see Figure 1). These two channels are separated by a small "monogroove" 
slot which creates a high capillary pressure difference and causes liquid to be 
pumped from the liquid channel to the circumferential grooves in the vapor 
channel. This configuration permits the axial transport and radial heat trans 
fer phases to be handled independently resulting in a high axial heat transport 
capability. 



The initial effort in the development of this heat pipe at NASA- JSC has 
been concentrated on a feasibility demonstration of che dual passage concept. 
Recently, investigation has been undertaken on the priming capabilities and 
behavioral characteristics of the liquid-vapor interface configuration during 
subjection to low-g or zero-g environments, similar to those which would be 
encountered during the operation of a low orbit Space Operation Center. 



Figure 1 Grumman Monogroove Heat Pipe 

This particular report formulates a mathematical model and computer pro- 
gram which describe the operational and performance characteristics of the 
Grumman dual passage heat pipe. This model allows us to predict the steady 
state performance when given the necessary physical parameters. This model 
can be used for the following: 

• To support the desigining and testing of laboratory test elements 
and prototypes. 

• To define the operating limits of the system. 

• To verify and correlate the data of element and prototype tests 
under 1-G and extrapolate them to 0-G environments. 

• To analyze the effect of varying heat source and sink temperatures 
on the thermal performance of the system. 

• To predict and simulate the thermal performance of a heat pipe 
radiator system operated under orbital environments. 





Computer modeling and simulation of heat pipes is a relatively new 
area with a majority of the work having been accomplished over the last 
10-12 years. Early modeling of heat p’pes was accomplished by S.W.Chi 
(2,3) with later contributions made by Thrush et al . (4) and D.K. Anand(5). 
Frank (6) has developed a generalized heat pipe equation and optimization 
method fcr grooved heat pipes which provides a methodology for determining 
optimum groove dimensions. Finally, Holm and Miller (7) completed a para- 
metric study of the defining equations for heat pipe operation to predict 
the performance characteristics from one which is dimensionally similar 
using similarity relations. 

The fundamental technique used in a majority of the works cited above 
was that of solving the pressure equations simultaneously, A similar approach 
has been adopted here. 


ANALYSIS 


Alario et al (1) have presented the two differential pressure balance 
relationships which govern the performance of a monogroove heat pipe. 
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Equation (A) assures that the wall wick capillary pressure rise is 
sufficient to overcome the total viscous pressure losses in the vapor channel, 
liquid channel and circumferential wall grooves, plus the gravity head losses 
associated with the inside diameter of the vapor channel along with any eleva- 
tion difference between the evaporator and condenser. In this particular 
investigation the primary concern lies in the performance during zero-g opera- 
tion resulting in a simplification of this equation. 

Equation (B) examines the pressure change resulting from the monogroove 
slot to insure that sufficient pressure is developed, to overcome the vapor 
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and liquid viscous losses, plus the gravity loss due to adverse tilt. Simul- 
taneous solution of equations (A) and (B) provides a method for determining 
the maximum heat transfer, Q. Equations (1) through (8) are as presented by 
Alario et. al. and describe the individual pressure difference terms as a 
function of the physical geometry. 
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Limitations on the heat pipe can be determined as follows: 


"When both of the differential pressure relationships are 
satisfied, the maximum heat transfer is governed by the 
wall wick capillary structure and the heat pipe performance 
is wall -wick limited. However, if the monogroove slot 
cannot sustain the necessary capillary pressure rise (i.e., 
large slot gap) then the heat transfer is prematurely 
limited." 

THE PROGRAM 

Essentially, the program is designed to solve equations (1) through (8) 
as a function of Q, substituting those values into equations (A) and (B) and 
then solving equations (A) and (B). 

Various checks are made throughout the program to determine if the heat 
pipe will prime properly in zero-g, if the flow is laminar or turbulent and 
if the sonic limit or entrainment limits have been exceeded. 

Table 1 shows the necessary input parameters, while Table 2 lists the 
computer nomenclature used in the program. A strong effort was made to use 
nomenclature consistent with that used by Grumman in their work in order to 
avoid confusion. This was accomplished by using a nomenclature listing ob- 
tained through Paul Marshal of NASA. In some instances it was necessary to 
deviate from this list either because the author was unaware of the existance 
of the Grumman term or a different approach was used requiring additional 
terms. 

Figure 2 illustrates a flow chart of the program and is self explanatory 
Numerous comments are included throughout the program which should help 
clarify the various steps, and an effort was made to structure the program in 
an orderly fashion. Appendix A contains a complete listing of the program, 
while Figure 3 is a copy of the output for a specific trial run. 

PROGRAM CAPABILITY 

Based upon the input physical parameters such as evaporator length, con- 
denser length, vapor passage diameter, liquid passage diameter etc., this pro 
gram is capable of determing the maximum heat transfer capacity in watts and 



and the transport capacity in watt-meters. In addition, the user can deter- 
mine if the sonic limit or entrainment limits have been reached, whether 
the heat pipe is wall wick or monogroove limited, if the heat pipe will prime 
properly in zero-g and the comparative values of the gravity head, monogroove 
head and the net capillary rise. 

Through the use of this program and a simple incremental loop,*** the 
effect of each of the input parameters can be determined individually answer- 
ing such questions as "What happens if the evaporator length is increased?" 
or "What if the vapor diameter is decreased?" This technique although not a 
sophisticated optimization technique, can provide valuable information as to 
the importance of the various physical parameters. 

PROGRAM VERIFICATION 

Preliminary testing of the program was accomplished by comparing computer 
predictions with actual experimental results available to the author. Figure 
3 is an example of one such trial run. The experimental results of actual pro- 
totype tests were very limited, but initial verification tests indicate that 
the program predicts the heat transfer capacity with reasonable accuracy with 
all deviation from actual results occurring on the low side. That is to say 
the program underestimates the actual capacity to a small degree. 

In order to determine an actual percentage error in the computer predic- 
tions, additional experimental data would need to be obtained from Grumman and 
comparisons made. 


*** This has been done and will be discussed in the oral presentation. 



TABLE 1 : INPUT VALUES 


PROPS (12) 

Fluid Properties at OPTEMP 


XLAT 

Latest Heat of Vaporization 

KJ/kg 

RHOL 

Liquid Density 

Kg/m^ 

RHOV 

Vapor Density 

kg/m^ 

TCONL 

Thermal Conductivity Liquid 

W/m°C 

XMUL 

Liquid Viscosity 

Centi Poi 

XMUV 

Vapor Viscosity 

Centi Poi 

PSAT 

Saturation Pressure 

Bar 

CP 

Specific Heat Constant P 

KJ/kg °C 

STEM 

Surface Tension 

N/M 

TCRIT 

Critical Temp 

°C 

PCRIT 

Critical Pressure 

Bar 

RMW 

Gas Constant X Molecular Weight 

kj/kg °K 


GEOMETRY 

EVAPM 

Evaporator Section Length 

M 

TRANSM 

Transfer Section Length 

M 

CONDM 

Condenser Section Length 

M 

DV 

Diameter Vapor Tube 

MM 

DL 

Diameter Liquid Tube 

MM 

TW 

Wall Thickness 

MM 

YS 

Yield Strength 


AZ 

Liquid Area Fraction 

0/ 

IQ 

THETWR 

Fluid Wetting Angle 

Radians 

TILT 

Tilt Height 

MM 

TW 

Wall Thickness 

MM 


TABLE 2: COMPUTER NOMENCLATURE 


GEOMETRY 


LEFFM 

Effective Length 

M 

LVM 

Overall Length 

M 

IXX 

Index 


WW 

Dummy Variables Web 


RW 

Root 


GD 

Groove 


XX 



SLPHTR 

Taper Angle of Wall Wick 

Radians 

AW 

Wetted Area Dummy 


WPW 

Wetted Penmeter Dummy 


DW 

Wetted Diameter Dummy 


AWE 

Wetted Area Evaporator 

m2 

WPWE 

Wetted Perimeter Evaporator 

M 

OWE 

Wetted Diameter Evaporator 

M 

AWC 

Wetted Area Condenser 

m2 

WPWC 

Wetted Perimeter Condenser 

M 

owe 

Wetted Diameter Condenser 

M 

OD 

Outside Diameter 

MM 

ODM 

Outside Diameter 

M 


HYDRAULIC 

DIAMETERS Change from to 

WPL 

Wetted Perimeter Liquid 

AL 

Area Liquid Channel 

DLH 

Diameter Liquid Channel 

WPV 

Wetted Perimeter Vapor Channel 

AV 

Area Vapor Channel 

DVH 

Diameter Vapor Channel 



VARIABLES 


WALL WICK 


WWE 

Web Width Evaporator 


RWE 

Root Width Evaporator 

m 

GDE 

Groove Depth Evaporator 

m m 

TRIE 

Thds/inch Evaporator 

thd/in 

NE 

Number of Evaporator Feeds 


WWC 

Web Width Condenser 

m m 

RWC 

Root Width Condenser 

m m 

GDC 

Groove Depth Condenser 

m m 

TPIC 

Thds/inch Condenser 

Thds/in 

NC 

Number of Condenser Feeds 


QZ 

Form Factor for Heat Transfer 



MONOGROOVE 


WM 

Monogroove Width 

m m 

ALPAD 

Monogroove Taper Angle 

Degrees 

DELT 

Monogroove Stand Off Distance 

m m 

GAMMA 

CP/ 

CV 

m/s^ 

G 

Acceleration of Gravity 


HEAD LOSS COEFFICIENTS 


RL 

fc/p.'iX Dummy Variable 


CA 

Liquid Channel LC 

ks o 

kj • ms 

CWE 

Wall Wick LC Evaporator 

II 

CWC 

Wall Wick LC Condenber 

II 

CW 

Wall Wick LC Total 

II 

PCW 

A P Wall Capillary 

N/m^ 

PCM 

A P Monogroove 

II 

PHEAD 

A P Head Dia.pt DVH 

II 

PTILT 

A P Tilt Tilt 

H 

XLMG 


-- 

CMGR 

Monogroove Radial LC 

-- 

CMG 

CMG ^ CMGR CMGR ^ 0.0 


CVL 

Laminar Vapor Channel LC 

kg o 
kj'm's" 

RV 

v/ pv# A Dummy Variable 


EQNA 



CTOT 

Rt Side of E 'NA 


PNET 

Lft Side of EQNA 

N/m^ 

Q 

Heat Transport 

KJ/s 

EqNB 



PCMN 

Lft Side of EQNB 

N/m^ 

PVL 

Rt Side of EQNA 

N/m^ 

REYNOLD"S NO. 



REN 1 

D/ XvA • A 


REN 

Reynold's Number REN 1*Q 

DQ 

CVT 

Turbulent Loss Coefficient 

. Kg 2 

QW 

Heat Transport 

ik»m*s 

Matts 

RCWM 

Transport Capacity 

W • m 

QSLW 

Heat Transport at *^onice Limit 

Watts 

TCSLW 

Transport Cap at Sonic Limit 

W • m 

QELW 

Heat Transport @ Entrainment Limit 

Watts 

TCELW 

Transport Cap @ Entrainment Limit 

W • m 



FIGURE 2 FLOW CHART 
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PRINT OUT LOSS COEFFICIENTS 
AND HEAT TRANSPORT RESULTS 
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C c 

C DUAL PASSAGE HEAT PIPE ANALYSIS 

C 

C HM 11/18/82 

C 

C 


DIHENSION PR0P<12) 

REAL LEFFH»LVH»NC»NE»OTEHPC 
PI = 3.1415926 

OPEN < UNI T=1.NAME= 'HEAT. OUT ' »TYPE=' UNKNOWN' > 


C 


C 

C 


T 

I 


i 


4' 


I 


I 


I 



WRITE (1»10) 

10 FORMAT (' '»/»20Xr 'MATHEMATICAL MODELLING OF HEAT TRANSFER IN 
+ DUAL PASSAGE HEAT PIPES') 

WRITE(lf20) 

20 FORMAT(' 25X »' **»*»»»#♦*#»#♦* OUTPUT DATA ♦♦****»*»»♦»»*♦»*'./) 


ID = 1 

CALL FLUID(PROP»OTEMPCf ID) 
XLAT=PR0P(1) 

RH0L=PR0P(2) 

RHOV=PROP(3) 

TC0NL = PRC)P(4) 

XMUL=PR0P<5) 

XMUV=PR0P(6) 

PSAT--PRQP(7) 

CP =PR0P(B) 

STEN=PR0P',7' 

TCRIT=PROP( 10) 

PCRIT = PROPn.l ) 

RMW= PROP a 2) 

DATA STEMPCr G / 40, Of 9.81 / 


PRINT OUT FLUID PROPERTIES 


URITEClf 100) 

100 FORMAT(' ' f20Xf '»*»** FLUID PROPERTIES *♦»♦)»:') 
WRITE(lfllO) OTEMPCfSTEMPC 

110 FORMAT(' '.'OPERATING TEMP = 'fT20fF10.4f C '.T40. 
+ 'STORAGE TEMP = ' . T60 . F 10 , 4 . ' C ) 

WRITE (lfl20) rCRITfPCRIT 

120 FORMAT<' '.'CRITICAL TEMP = ' .T20.F10.4. ' C .140” 

+ 'CRITICAL PRESSURE = ' » T60.F10.4. ' C ') 
WRITE(1.130) XLAT.TCONL 

130 FORMAT<' '.'LATENT HEAT = ' . T20 . FIO . 4 . ' KJ/KG'.T40. 
+ 'THERMAL COND LIQ = ' . T60 » FI 0 . 4 . ' W/M*C ' ) 


WRITE(1.140) RHOL.RHOU 

140 FORMAT(' '.'LIQ DENSITY = ' . T20. FIO . 4 . ' KG/M”3 '.T40. 

+ 'VAP. DENSITY = ' . T60 . F 10 . 4 . ' KG/M''3') 

WRITEd.150) XMUL.XMUV 

150 FORMAT(' '.'LIQ. VISCOSITY = ' . T20 . FIO . 4 . ' CPOISE '.T40. 
+ 'VAP. VISCOSITY = ' .T60.F10.4. ' CPOISE') 

WRITE(1.160) CP.STEN 

160 FORMAT(' '.'VAP SPEC HEAT = ' . T20.F10 . 4. ' KJ/KG*C '.T40. 
+ 'LIQ SURF. TENSION = ' . T60 .FIO . 4. ' N/M ') 

WRITE(1.170) PSAT.RMW 

170 FORMAT(' '.'SAT. PRESSURE = ' . T20 . FIO . 4 . ' BAR '.T40. 

+ 'SPEC. GAS CONST = ' . T60 . FIO . 4 . ' KJ/KG*K ') 
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INPUT GEOMETRY 
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DATA EVAPMf TRANSM»CONDM»DV»DLf TU /O . IS » 0 . 46 f 0 * 15 » 13 . 39 » 6 . il' » 1 . 24/ 

C YIELD STRENGTHf LIQUID AREA FRACTIONr WETTING ANGLf 
DATA YS»AZ»THETWR»TILT / 1 37 . 9E+06 » 1 » 0 . 0 . 3 . 1 75/ 

C 

C INPUT WALL WICK PROPERTIES 

C WEB WIDTHf ROOT WIDTHf GROOVE DEPTH» THREADS/INCH, NO. OP FEEDS | 

DATA WWE»RWEfGDE»TPIE»NE /O ♦ 051 » 0 . 014 »0 ♦ 196 » 192. Of 1 . 0/ # 

DATA WWCfRUCfGDCf TPlCfNC /0 . 051 f 0 . 014 f 0 . 1 96 f 1 92 . 0 1 1 . 0/ 

C ;-ORM FACTOR FOR HEAT INPUT 
DATA QZ /l.O/ 

C 

C INPUT MONOGROOVE PROPERTIES 
C GROOVE WIDTHf TAPERf STAND OFF DISTANCE 
DATA WMf ALPHADfDELT /O . 254 f 0 . 0 f 1 . 24/ 

CONVERT INPUT VARIABLES FROM (MM) lO <M) 

DV = DV/1000. 

DL = DL/1000. 

TW = TW/1000. 

WWE = WWE/1000. 

RWE = RWE/1000. 

GDE = GDE/1000. 

WWC = WWC/1000. 

RWC == RWC/1000. 

GDC = GDC/ 1000. 

WM = WM/1000. 

DELT = DEL f/ 1000. 

TILT = Til. r/1000, 

CALCULATE 6E0METf<Y 

1. PIPE GEOMETRY 
LEFFM = TRANSM+(EUAPM + L'ONDM> /2 
LVM = TRANSM + EVAPM + CONIii 

2. WALL GROOVE GEOMETRY (S.'nPUFlED) 

FOR A FIRST APPROX. ASSUME GROOVE FILLED WITH LIQUID 
(IE. AWSEG=0.) 

IXX = 0 
Wl* = WWE 
RW = RWE 
GD = GDE 

40 XX = (WW-RW)/(2*GD) 

ALPHTR = ATAN(XX) 

AW = GD*(WW+RW)/2.0 
WPW = RW+(2*GD) /COS (ALPHTR) 

DW = (4«AW)/WPW 
IF (IXX .NE. 0) GO TO 50 
AWE = AW 
WPWE = WPW 
DWE = DW 
WW = WWC 
RW = RWC 
GD = GDC 
IXX = 1 
GO TO 40 
50 AWC = AW 
WPWC = WPW 
DWC = DW 

3. monogroove; geometry 

OD = DV + DL + (2*TW) + DELT 

4. HYDRAULIC DIAMETERS 

ASqfJMP A7 K 1 { U\nv FIIM ^ 
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C A. LIQUID HYD DIA 

UPL = PI*DL»A2 
AL = WPL»DL/4 
DLH = DL 

B. VAPOR HYD* DIA. 

UNTIL DEFINITION OF NE IS KNOWN 
WPV = PI *00 

AV = WPV»DV/4 _ 

OVH = DV ORIGINAL PAGE IS 

OF POOR QUALflY 

CHECK IF DVH > DLH 

IF <DVH.GT.DLH) GO TO 70 
WRITE(lf60) 

TYPE 60 

60 F0RMAT<' ./,'***»♦ ERROR-DIAMETER OF LIQUID CHANNEL EXCEEDS 
+ DIAMETER OF VAPOR CHANNEL **♦*♦',/) 

70 CONTINUE 

CHECK PRIMING CAPABILITY OF VOLUME 

ALl = <PI*DL#DL/4.0)+<DLLT*UM) 

DLl = SQRT(4»AL1/PI ) 

AVI = PI»DV*liV/4.0 
AV2 = AVl-ALl 
DV2 = SQRT(4.0*AV2/PI ) 

IF (DLl .GT. DV2) GO TO 200 
WRITE (1»80> 

80 FORMAT (' ' » 20X » ' »***» WILL NOT PRIME JN ZERO-6 t**t*'> 

GO TO 6000 
RETURN 

200 CONTINUE 

TPME = TPIE*39.37 
TPMC = TPIC » 39.37 
CONVERT TPI TO THDS/M (TPM) 

PRINT OUT GEOMETRY 

OUTPUT HEAT PIPE GEOMETRY 
WRITE(1»300) 

300 FORMAT(' ' . // » 20X » ' HEAT PIPE GEOMETRY #*»*♦',/) 

WRITE(1»310) EVAPMf CONDM tTRANSM 

310 FORMAT(' '»'EVAP LENGTH = '»F10.4f' M COND LENGTH = '»F10.4, 

+' M TRANSFER LEN. = 'fF10.4f' M ') 

WRITE(1»320) DVH»DLH»OD 

320 FORMAT(' 'f'VAP CHAN HYD DIA == »F10.6»' M LIQ CHAN HYD DIA - 

+ F10.6»' M O.D. = fF10.6»' M ') 

WRITE(1.325) AV,AL 

325 FORMAT(' '»'AREA VAP CHAN ^ ',E14.6»' ^M”2) AREA LIQ CHAN = '» 

+ E14.6t' <M”2) '»/) 

WRITE( 1f330) WWE » TPME f WWC » TPMC 

330 FORMAT(' 'f'WALL WICK OPENING: EVAP = 'fF10.6f' Mf 'fF10.4f 
+ 'TPM COND = 'fF10.6f' Mf 'fF10.4f' TPM ') 

WRITE<1f340) UMfDELT 

340 FORMAT(' ' f ' MONOGROOVE OPENING = 'fF10.6f' Mf LENGTH = 'fF10.6f 
+ ' M ' f/) 


CALCULATE LOSS COEFFICIENTS 

1. CA - LIQUID CHANNEL LOSS COEFF. (LC) 

RL = XMUL/(RHOL»XLATH[1000, ) 

CA = <32#RL*LEFFM)/<AL*DLH*DLH) 

2. CW - WALL WICK LC 

. 3,2JCJ='I*RL*DVH*QZ 
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CUE = CW/(8*NE#EyAPH*TPME*AWE*DWE)(cDWE) 
cue = CU/<8HtNC*C0NDM)«TPME»AUC4HiUC»DWC) 
CW = CUE + cue 

3. PCU - DELTA P UALL CAPILLARY 

PCW = 2»STeN*C0S(THETWR+ALPHTR)/WUE 


4. FCM - DELTA P MONOGROOVE. 

PCM = 2»STEN*C03(THETWR)/UM 

5. PHEAD - DELTA P HEAD DIA. IN GRAVITY 
PHEAD = RHOL*DVH*G 

6. PTILT - DELTA P TILT IN GRAVITY 

PTILI = RHOL*TILT*G ORI3irOi's. - -*2 S!| 

CHECK FOR 0-G. (TILT < -100) OF POOR Q’^ALfTY 

IF (TILT .GT. -100.) GO TO 1000 
PHEAD = 0.0000 
PTILT = 0.0000 
1000 CONTINUE 

l HD = PHEAD + PTILT 


EONS 1»2»3»4»7 i»8 HAVE BEEN COMPUTED 

NOU IT MUST BE DETERMINED IF FLOW IS LAMINAR 

(FOR EQN 5) OR TURBULENT (FOR EQN 6) 

FIRST ASSUME LAMINAR FLOW. CALCULATE Q 
THEN CALC. REN AND CHECK. 

8. eVL - LAMINAR VAPOR CHANNEL LOSS COEFFICIENT 

RV = XMUV/(RH0V*XLAT»1000.0) 
eVL = (32»RV*LEFFM)/( AV*DVH*DVH> 

heat transport calculations -- 


COMPUTE Q FROM EQN A 
CTOT = CVL+CATCU 
PNET = PCU-PHD 
Q = PNET/CTOT 
NOW CHECK EQN B 
PCMN = PCM-PTILI 
PVL = (CA+CVL)*Q 
IF (PCMN .GE. PVL) GO TO 2000 
Q = PCMN/(CA+CVL) 

PNET - PCMN 

CALCULATE REYNOLD'S NUMBER 

2000 RENl = (DVH*i000.0)/ (XMUViKAViKXLAT) 

2300 REN = REN1*Q 
TYPE *f REN 
C 

C CHECK REYNOLD'S NO. FOR TURBULENT OR TRANSITION FLOW 
2500 IF (REN .LT. 2300.) GO TO 4050 
IF (REN .GT. 3900.) GO TO 3900 
C 

C TRANSITION FLOW CALCULATIONS 

FTRAN = (3.31E-5)»(REN»«.69) 

GO TO 4000 
C 

C TURBULENT FLOW CALCULATIONS 
3900 FTURB * .0791/(REN)K»,2500) 

FTRAN * FTURB 
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C TURBULENT VAPOR LOSS COEFFICIENT 
4000 CVT = (2*FTRAN)KLEFFM)/( AV*AV*UVH*RHOV*XLAT*XLAT) 


COMPUTE Q FROM EON A FOR TURBULENT FLOW 
IF (PNET .GT, PCMN) GO TO 4030 
B = CA/CVT 
GO TO 4040 
0 B=(CA+CW)/CVT 
0 XX = (BiKB)+4*(PNET/CVT) 

XX = SQRT(XX) 

Q1 = (XX-B)/2.0 

IF DIFFERENCE BETWEEN NEW 0 (CAUSED BY TURB. FLOW 
CALCULATIONS) AND PREVIOUSLY CALCULATED 0 IS >1X 
THEN RECALC. REYNOLD'S NO. USING NEW Q AND REPFAI 
UNTIL DIFF < IX 


DIFF = ABS( (Ql-Q)/Q) 

Q = Q1 

IF (IIFF.GE. .010) GO TO 2300 
CHECK EON B 

PVLT = Q*(CA+(CVTiKQ) ) ORIGINA2- PAGE 13 

IF (PCMN ,GE. PVLT) GO TO 4050 Op POOR QUALFIY 

ZY = (CA»CA) + 4.0*CVTHtPCMN 
ZY = SORT(ZY) 

0 = (ZY-CA)/(2.0»CVT) 

PVu = PVLT 
PNET = PCMN 
GO TO 2300 
4050 CONTINUE 


REPEAT LAMINAR EQUATIONS 
REN = RENDKQ 

IF (REN .BT. 2300.) GO TO 4060 
PVL = (CA+CVL)*Q 
IF (PCMN .6E. PVL.) GO TO 4060 
0 = PCMN/(CA+CVL) 

PNET = PCMN 
4060 CONTINUE 

QW = Om 1000.0 
TCWM = QWHtLEFFM 

END HEAT TRANSPORT CALCULATIONS 


CHECK SONIC LIMIT 


GAMMA = 1.33 

TEMPK = OTEMPC -F 273 

VS = SQRT(1000.0#GAMMA»RMW)KTEMPK/(2.0)*;GAMMA + 2.0) ) 

QSLW = RHOV*XLAT*AV*VS»1000.0 
SLTC = QSL * LEFFM 

CHECK ENTRAINMENT LIMIT 

XX = (RHOV*STENfXLAT*XLAT)/(WM/(2*COS(THETWR) ) ) 

QELW = AV#SQRT(XX)*1000.0 
TCELW=QELW*LEFFM 

PRESSURE AND LOSS COEFFICIENT OUTPUT 

WRITE (lf400) 

400 FORMAT(' ' » // » 20X » ' **»*» LOSS COEFFICIENTS ( KG/K J*M*S''2 ) t 

WRITE (1»410) 

410 FORMAT(' '» 'LIQUID CHANNEL. VAPOR CHANNEL WALL WICK 

■F WALL WICK EVAP WALL WICK COND') 

IF CREN.QF.2300. ) fiO TO 5000 . . _ ______ 



WRITE (lf420) CAfCVLfCWfCWEfCWC 
GO TO 5100 

5000 WRITE <1»420) CA*CVr»CW»CWE.CWC 
5100 CONTINUE 

420 FORMAK' ' » El 4 . 6 » 4X » E14 . 6 » 2X » El 4 . 6 » 2X » El 4 . 6 » 4X i El 4 . 6 ) 

WRITE <1»430) 

430 FORMAT(' ' » / » 20X » ' ♦»»*» CAPILLARY PRESSURES < KU/MikS'':.' ) ****'f/) 
WRITE <1»440) 

440 FORMAT (' 'r'WALL GROOVES MONOGROOVE TOTAL GRAVITY HEAD 

+ NET CAP RISE NET MONOGROOVE') 

WRITE (1»450) PCW»PCM»PHD»PNET»PCMN 
450 FORMAT(' ' » FI 0 . 4 f 5X » F 1 0 . 4 » 8X » F 10 , 4 » 1 2X » F 10 . 4 » 9X » F 10 ♦ 4 f // ) 

WRITE <1»460) REN 

460 FORMAT (' '» 'REYNOLDS NUMBER = 'fFlO.4) 

C 

c 

C PRINT TRANSPORT CAPACITY 

C 

WRITE (1»490) TCWM»QW 

490 FORMAT (' 'TRANSPORT CAPACITY '»F10.2»' WATT#M OMAX = 

+ F10.2f ' WATTS ' ) 

IF (PCMN.GT.PVL) GO TO 5200 
WRITE (1»500) 

500 FORMAT! '+'»' ***** MONOGROOVE LIMITED *****') 

TYPE *f '***** MONOGROOVE LIMITED *****' 

GO TO 5300 original page wi 

5200 WRITE (1»510) nF POOR QUALfTY 

510 FORMAT('+'»' ***** WALL WICK LIMITED ***** , 

TYPE *f '***** WALL WICK LIMITED *****' 

5300 IF(QW .LE. QSLW) GO TO 5400 
WRITE <lf520) 

520 FORMAT! ' + '»' He**** SONIC LIMIT EXCEEDED ***** ') 

WRITE !1»530) QSLW»SLTC 

530 FORMAT!' ' » 20X » ' L IMI TS ARE: Q == ' » T30 » F ] 0 . 2 r ' WATTS'»T60» 

+ 'TRANSPORT CAPACITY = ' f T80 f F 1 0 , 2 r ' WATT-M') 

5400 IF !QW .LE, QELW) 60 TO 6000 


DATA FORMAT 


WRITE !1»540) 

540 FORMAT! ***** ENTRAINMENT LIMIT EXCEEDED ***** '; 

WRITE !lf530) QELWfTCELW 
6000 CONTINUE 

CLOSE !UNIT=1) 

STOP 

END 

C 7 

C FLUID PROPERTIES SUBROUTINE 

C 

C 

c 
c 
c 

C *I 

c * 

C *F 

C # 

C * 

C * 

C * 

c * 

c * 

c * 

c * 

c « 

c « 

c * 

t 


!l-10f 

IfX) - TCRIT f 

PCRIT.TMPi 

TEMP 

-60»-40f-20»0 

»20.40»60i 

IPERTY 

DAT!X»#»X) 

PROP' 

XLAT 

2 

1 

RHOL 

3 

2 

RHOV 

4 

3 

TCONL 

5 

4 

XMUL 

6 

5 

XMUV 

7 

6 

PSAT 

8 

7 

CP 

9 

8 

STEN 

10 

9 

TCRIT 

1 rl 

10 

PCRIT 

2»1 

11 

1 

3.1 

1 

'rr Vi /-. rti 


UNITS 

KJ/KG 

KG/M3 

KG/M3 

KJ/M DEG.C 
CENTIPOISE 
CENTIPOISE 
BAR 

KJ/KG DEG.C 
N/M 

DEG. C 
BAR 

juaa-'f'p-'’ • 
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FLUID*'? * ETC. 

BLOCK DATA 

DIMENSION AhMCNClOflO) »DAT(10»10»1 ) 

DIMENSION PRCPSl (10) fXLATl ( 10) f RHOLl ( 10) fRHOVl ( 10 y »TC0N1 (10) 
DIMENSION XMULl ( 1 0 ) j XMUV l ( 10 ) » PSAT 1 ( 1 0 ) » CPI ( 10) >STEN1 (10) 

COMMON /AMMON/ PROPSl > XLATl- RHOLlf RHOOlf 
+ TCONlf XMULl » XMUV1» PSATlr CPI » STEM 
EQUIVALENCE ( AMM0N( 1 f 1 ) f PROPSl ) 

EQUIVALENCE ( DAT ( 1 . 1 > 1 )» AMMON ) 

ENTER DATA IN THE FOLLOWING BLOCK ACCORDING TO ABOVE FORMAT 

2345 / STATEMENT FIELD (72 SPACES) 

AMMONIA PROPERTIES 


DATA PROPSl /132.4»112.9.-77,7»--33.4» .4882»-60»120»0f0,0/f 
XLATl /1434. f 1384. » .338 . » 1263. » 1187. f 1101 . »1026. r 891 . . 699 . . 428 . / » 


RHOLl 


1.505.7.455.1 . 


RHOVl 

TCONl 


/71 4. 4. 690. 4. 665. 5. 638. 6. 61 0.3. 579. 5. 545.: 

374.4/. 

/. 03. .05. 1.62.3.48.6.69. 12. 0.20. 49 i 34. 13.54.92. 113, 16/. 

/. 294. .303. .304. .298. .286. .272. .255. .235. .212. .184/ 
DATA XMULl / . 360 . . 290 . . 260 . . 250 . . 220 . . 200 » . 1 70 . . 1 50 . . 1 1 0 . , 70/ . 
XMUVl /.0072. .0079. .0085. .0092. , 0 1 0 1 . . 00 16 . . 0 1 27 . . 0 1 4 • . 0 1 6 . 
.0189/. 

PSATl /. 27. .76.1.93.4.24.8,46.15.34.29.80. 40.9. 63.12. 90.44/. 
CPI /2. 05. 2. 075. 2. 1.2. 125. 2. 15. 2.16. 2.18. 2.21 . 2.26. 2.92/ 
DATA STENl /. 04062. .03574. .0309. .0248. . 021 33 . . 0 1 833 . . 0 1 367 . 
.00767. .005. .0015 / 


ACETONE PROPERTIES 


END 


FLUID PROPERTY CALCULATION SUBROUTINE 


SUBROUTINE FLUID ( PROP . OTEMPC . I D ) 

DIMENSION DAKlO.lO.l) . PR0P(9). AMMON(IO.IO) 

COMMON /AMMON/ XI ( 10 ) . X2 ( 10 ) . X3 ( 10 ) . X4 ( 10 ) . X5 ( 1 0 ) . 

+ X6(10) .X7(10) .X8(10) .X9(10) .XIO(IO) 

EQUIVALENCE (DAT. AMMON) 

EQUIVALENCE (AMMON. XI) 

C 

100 TYPE INPUT OPERATING TEMP IN DEG. C. ' 

120 ACCEPT ♦. OTEMPC 

TYPE ».'OTEMPC= '.OTEMPC 

TYPE 10. ( (DAT( I . J. 1 ) . 1=1 . 10) . J=1 . 10) 

10 FORMAT( ' ' . 10F8.3) 

C 

C CHECK TEMP TO SEE OTEMPC > TCRIT 

IF (OTEMPC .LT. DAT(l.l.ID)) GO TO 210 
TYPE ».'»**» WORKING TEMP EXCEEDS TCRIT **** ' 

GO TO no 
C 

C CHECK FOR TEMP > FREEZING 

210 IF (OTEMPC .GT. DAT(3.1.ID)) GO TO 250 

TYPE TEMPERATURE BELOW FREEZING POINT 

GO TO no 


r 





r 

I 


I 


i 

% 


i 

4 



I 


1 


C CHECK FOR OTEMPC IN RANGE OF DATA - TDMIN .. OTEMF'C . TDHAX 

250 IF (OTEMPC .LT. DAT(6»1*ID)) GO TO 200 
IF (OTEMPC .GT. DAT(7»1»ID)) GO TO 300 
GO TO 350 

300 TYPE tf'tt** INSUFFICIENT DATA - TEMP OUT OF DATA RANGE **■*$' 
TYPE *f DAT(6» I » ID) f ' - ' » D AT ( 7 r 1 f I D ) 

310 60 TO 110 
C 

C CALCULATE INDEX NO. 

350 XNO = (OTEMPC/20.) + 4. 

NMl = INT(XNO) 

N = NMl + 1 
R = N - XNO 
DO 400 I = lx? 

II I + 1 

PROP ( I ) DAT ( N f 1 1 . .1 D ) - ( R» (DAT M , 1 1 , 1 D ) -DAT ( NM 1 » 1 1 » I n ) ) > 

400 CONTINUE 

PROP (10) = DAT(1»1?ID) 

PROP(ll) = DAT(2flfID) 

PR0P(12) = DAKSjIjID) 

GO TO 450 

no TYPE 'INPUT NEW TEMP (DEG, C) ' 

GO TO 120 
450 CONTINUE 
RETURN 
END 


'} 

i 
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